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The linear O(N) sigma model undergoes a symmetry restoring phase transition at finite 
temperature. We show that the nonlinear O(N) sigma model also undergoes a symmetry 
restoring phase transition; the critical temperatures are the same when the linear model 
is treated in mean field approximation and the nonlinear model is treated to leading plus 
subleading order in the 1/N expansion. We also carefully define and study the behavior 
of f n and the scalar condensate at low temperatures in both models, showing that they 
are independent of field redefinition. 
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1 Introduction 



The 0(N) model as a quantum field theory in d+1 dimensions JT| is a basis or prototype 
for many interesting physical systems. The bosonic field $ has N components. When 
the Lagrangian is such that the vacuum state exhibits spontaneous symmetry breaking 
it is known as a sigma model. This is the case of interest to us here. In d = 3 space 
dimensions the linear sigma model has the potential 

i(* 2 -/- 2 ) 2 

where A is a positive coupling constant and f w is the pion decay constant. The model 
is renormalizable. In the limit that A — > oo the potential goes over to a delta function 
constraint on the length of the field vector and is then known as a nonlinear sigma model. 

The classical limit of the field theory is obtained by neglecting or freezing out the time 
variable, leaving a field theory in d dimensions. In this limit only the zero Matsubara 
frequency of the full d+1 dimensional theory contributes to the partition function, and 
the temperature acts like a coupling constant. One then has a description of an 0(N) 
Heisenberg magnet in d dimensions which is a model for real material systems. This 
subject has a vast literature [g, [J. 

When N = 4 one has a model for the low energy dynamics of Quantum Chromody- 
namics (QCD). More explicitly, it is essentially the unique description of the dynamics of 
very soft pions. This is basically due to the isomorphism between the groups 0(4) and 
SU(2) xSU(2), the latter being the appropriate group for two flavors of massless quarks 
in QCD. The linear sigma model, including the nucleon, goes back to the work of Gell- 
Mann and Levy ||. This subject also has a vast literature. In the last decade much 
work has been done on chiral perturbation theory which starts with the nonlinear sigma 
model and adds higher order, nonrenormalizable, terms to the Lagrangian, ordered by 
the dimensionality of the coefficients or field derivatives |J. This whole program really 
has its origins in the classic works of Weinberg || |7|. 

Finally, the standard model of the electroweak interactions, due to Weinberg, Salam, 
and Glashow has an SU(2) doublet scalar Higgs field responsible for spontaneous sym- 
metry breaking. If one neglects spin-1 gauge fields the Higgs sector is also an 0(4) field 
theory. 

All of these limits are interesting to study at finite temperature. Magnetic materials 
typically undergo a phase transition from an ordered to a disordered state. If quarks 
are massless, QCD is expected to undergo a chiral symmetry restoring phase transition 
H, EJ . This may have implications for high energy nucleus- nucleus collisions; see especially 



TO, 11] in this respect. The electroweak theory is expected to have a symmetry restoring 



phase transition, too, at which point the baryon number of the early universe would have 
been finally determined [ |i~2|j . 

The linear sigma model was studied in the classic papers on relativistic quantum field 
theories at finite temperature [13], O, |T5||. The usefulness of leaving N as a parameter 
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arises from the fact that there is only one other parameter in the problem, the quartic 
coupling constant A (f n just sets the scale and is held fixed in our considerations). For 
QCD at least, and perhaps for electroweak theory too (this is not known, it is related to 
the Higgs mass) the appropriate limit seems to be A 3> 1, possibly even infinity. This limit 
is the nonlinear sigma model. The only proposed expansion parameter for this model is 
1/N. For N = 4 the first few terms in this expansion may not be quantitatively reliable 
but it is a good start. At least in the theoretical world we can imagine N as large as 
we wish. Presumably the physics does not change qualitatively with N as long as it is 
greater than one. 

Our basic physical interest in this paper is QCD. Among the questions that are rou- 
tinely asked are: Does QCD have a finite temperature phase transition? If so, is it 
associated with color deconfinement, or with the restoration of the spontaneously broken 
chiral symmetry, or are they inextricably intertwined? What would be the order of this 
phase transition? What would be its critical temperature? How do pions, the Goldstone 
bosons of QCD, decouple as the temperature is raised? How do the quark and gluon 
condensates behave as functions of temperature? Restricting our attention to Nf flavors 
of massless quarks, there are strong arguments H, 0] and numerical lattice computations 



Tqj which say that for Nf = 2 there is a second order phase transition and for Nf > 3 



there is a first order phase transition. Lattice calculations also suggest strongly that chiral 
symmetry is restored and color is deconfined at the same temperature |17| . 

The QCD Lagrangian is invariant under XJ(Nf) xXJ(Nf) transformations. This is iso- 
morphic to SUi(A r j)xSU^(A^j)xUbaryon ( 1) x Uaxiai (1) ; that is, left- and right-handed chi- 
rality transformations of the quark fields, baryon number conservation, and the famous 
axial U(l) symmetry. The axial U(l) is broken by quantum effects, particularly instan- 
tons. Baryons can be added to the sigma model if desired, but we shall not do so. The 
O(N) sigma model has no vestige of the axial U(l), although there do exist versions 
of the sigma model based on other groups in which it can be incorporated. In general 
SU(iV/)xSU(iV/) is isomorphic to 0(iV|) only for Nf = 2. Otherwise the interactions 
are different. This restricts any potential quantitative results of our analysis of the O(N) 
model to QCD to iV = 4. The very interesting issue of whether the axial U(l) symmetry 
is restored at high temperatures or not cannot be addressed here 



One must be careful in understanding how sigma models are being applied to the 
study of QCD at finite temperature. At very low temperatures one may argue that the 
only degrees of freedom which are excited are pions. One may then use chiral perturbation 
theory to study the thermodynamic properties, as in the classic works of Leutwyler [TJ5[ . 
At the very lowest energies this is just the nonlinear sigma model. Note that in this 
domain one is really studying a full 3+1 dimensional quantum field theory. On the other 
hand, near the critical temperature one may argue || [| that the soft long wavelength 
modes of QCD are in the same universality class as the 0(4) Heisenberg magnet in 3 
spatial directions. Now one is studying a classical field theory. The parameters of the 
effective free energy functional near the critical temperature may not be simply related to 
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the parameters of the sigma field theory model at very low temperature. Support for this 
point of view comes from lattice QCD computations which show that critical exponents 
are consistent with those of the 0(4) Heisenberg magnet in 3 spatial dimensions |20f . 
Conversely, it has been pointed out that in certain models with composite bosons these 
arguments might break down [^1] because the compositeness is important near a second 
order phase transition. A particularly useful observation is that long range correlation 
functions may be dominated by the soft modes (pions, sigma meson, ...) near the critical 
temperature but the equation of state itself is not |TB[ ; it is dominated by the myriad 
of other degrees of freedom (Hagedorn or Particle Data Book compendium of mesons and 
baryons or all of the colored quarks and gluons). 

We would like to shed some light on just a few of the issues relating to the above 
discussion. First, from the outset we restrict our attention to the linear and nonlinear 
sigma models in 3+1 dimensions. Within these confines we address three specific topics: 
the existence and nature of a chiral symmetry restoring phase transition, the low temper- 
ature dependence of the "pion decay constant", and the low temperature dependence of 
the "quark condensate". The first of these involves folklore, and the answer, at least in 
the nonlinear sigma model, is either obvious to the reader or else very surprising. For the 
latter two we point out some popular misconceptions and reproduce the existing results 
(when N = A) while showing that they are invariant under field redefinition. 

First consider the linear sigma model. At zero temperature the effective potential has 
a shape similar to the bottom of a wine bottle. It is minimized by a nonzero value of the 
field; this is the condensate. As the temperature is increased, the radius of bottom of the 



potential shrinks, and goes to zero at a critical temperature of T c = y 12/ (N + 2) It is 
a second order symmetry restoring phase transition. In the nonlinear sigma model |$| 2 is 
fixed at the value f%. Therefore, it would seem, chiral symmetry breaking is built into the 
Lagrangian and there is no possibility of restoring it at finite temperature. Another way of 
saying this is that there is no order parameter which can go to zero at finite temperature. 
At least this is the folklore in much of the nuclear and particle physics community! On 
the other hand, the critical temperature in the linear model is independent of A in the 
mean field approximation, so one can take the limit A — > oo and still have a phase 
transition! The counter argument to this is that the phase transition can go away in the 
limit and so nothing special happens at the aforementioned value of T c in the nonlinear 
model. We shall study the nonlinear model directly in section 2 at finite temperature in 
the large N approximation. To leading order in N we shall show quite straightforwardly 
that despite the constraint the nonlinear model has a second order phase transition at a 
critical temperature equal to that of the linear model. We shall show that this persists to 
the next to leading order in N, although here we must make an additional high energy 
approximation. The order parameter is identified as is the nature of the two phases. 
Next we study what is meant by "the pion decay constant on finite temperature" . At 
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zero temperature a common definition is 



(o\ai\Ap)) = iUP^ 



ab 



which relates it to the matrix element of the axial vector current of QCD between the 
vacuum state and a one pion state of momentum p. It is difficult, though perhaps not 
impossible, to generalize this definition to finite temperature. Within the linear sigma 
model one sometimes sees in the literature fw(T) identified with the thermal average of the 
sigma field, v = (a) = | (*&)[, which is the radius of the bottom of the effective potential. 
In the nonlinear sigma model this radius is necessarily fixed at f n . How then can one 
understand the result of Gasser and Leutwyler p3[ 



r 2 



12 f I 



+ 



(2) 



which was obtained at low temperature in chiral perturbation theory? (At this order 
chiral perturbation theory and the nonlinear sigma model are the same.) This issue is 
addressed carefully in section 3. 

Finally, a quantity of much interest, especially for the application of QCD sum rules 
at finite temperature, is the temperature dependence of the quark condensate. Gasser 
and Leutwyler p3[ and also Gerber and Leutwyler [7M computed this quantity at low 



temperature to be 



(qq) = (0\qq\0) 



1 / T 2 



N 



f 



Nf-l 



T 2 



12/= 



+ 



(3) 



This is obviously a different temperature dependence than that of f n {T). In addition, 
one wonders how it is possible to obtain information on quark condensates from a theory 
which has no explicit reference to quarks? These and related topics are studied in section 
4. 

Before beginning the technical part we remark that we will deal with vacuum di- 
vergences in a simple way: we ignore them. That is, we will always drop temperature 
independent divergent quantities from our expressions. The linear sigma model is renor- 
malizable and this neglect can easily be rectified; it does not change any of the principle 
results. The nonlinear model is not renormalizable, although one may consider it to be 
the A —>■ oo limit of the linear model. The short distance physics of these models in the 
context of QCD is not correct anyway. 

Put another way, it is known that the partition function can be expressed in terms of 
the vacuum scattering amplitudes, or S-matrix elements, for arbitrary reactions involv- 
ing n particles going in and m particles coming out |25]. This is the relativistic virial 
expansion. Consider, for example, a real scalar field with a quartic interaction [26]. The 



one loop contribution to the partition function is just the free Bose gas expression. The 
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two loop contribution corresponds to two particle scattering with the amplitude evaluated 
at the tree level. The three loop contribution corresponds to two particles in and four 
particles out, plus three particles in and three particles out, plus four particles in and 
two particles out, with all scattering amplitudes evaluated at tree level. The three loop 
contribution also has a part which corresponds to a vacuum one loop correction to the 
two particle scattering amplitude. By dropping all vacuum loop divergencies we are doing 
a virial expansion with the S matrix evaluated at the tree level. 

An early systematic study of various approximation schemes for the linear sigma 
model at finite temperature is ^7j. An overall introduction to relativistic quantum field 
theory at finite temperature is [28|. A somewhat analogous study of the lattice O(N) 
Heisenberg model is [B9[ 



2 Chiral Phase Transition 

It is well known that the linear O(N) sigma model in 3+1 dimensions has a second order 
phase transition when treated in the mean field approximation. We briefly repeat that 
analysis here as a warmup and precursor to the study of the nonlinear sigma model which 
is not nearly as well studied. We must emphasize that the direct quantitative applicability 
to QCD is limited by at least two factors. The first is: These sigma models do not have 
quark and gluon degrees of freedom and so one can never describe high temperature 
quark-gluon plasma with them. In addition, the compositeness of the bosons may even 



influence the phase transition itself The second is: The group SU(iV/)xSU(iV/) is 
isomorphic to O(Nj) only for Nf = 2 and this limits the analogy to two flavors of massless 
quarks. Indeed, lattice computations of the nonlinear sigma model in 3 dimensions based 



on other groups show first order behavior [30 



The conventions and notation used here are consistent with those of 28 . 



2.1 Linear sigma model 

The linear sigma model Lagrangian is 

C = \{d^f-\(&-fl) 2 (4) 

where A is a positive coupling constant. The bosonic field $ has iV components. Rather 
arbitarily we define the first N — 1 components to represent a pion field 7r and the last, 
iV'th, component to represent the sigma field. Since the O(N) symmetry is broken to an 
0(N — 1) symmetry at low temperatures, we immediately allow for a sigma condensate 
v whose value is temperature-dependent and yet to be determined. We write 

$i(x,t) = 7T 4 (x,t) i = l,...N-l 

$jv(x,t) = v + a{x,t). (5) 
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In terms of these fields the Lagrangian is 

C = ~ (^tt) 2 + \ {d,af - \ (, 2 - fl + 2va + a 2 + tt 2 ) 2 . (6) 

The action at finite temperature is obtained by rotating to imaginary time, r = it, and 
integrating r from to (3 = 1/T. (However, we keep the Minkowski metric; hence, 
= djdx^ with 8q = d/dt = id/ dr.) The action is defined as 



S = -±(fi- v ff3V + J Q drJ v d"x[\[{d^f-^ + {d,af-mla 

+ \v (, 2 - fl) a - \va (tt 2 + a 2 ) - \ (a 2 + tt 2 ) 2 } , (7) 



where the effective masses are 



K = A(V-/ 2 ) 
K = A(3^ 2 -/ 2 



At any temperature v is chosen such that (a) = 0. This eliminates any one particle 
reducible (1PR) diagrams in perturbation theory, leaving only one particle irreducible 
(1PI) diagrams. 

At zero temperature the potential is minimized when v = f n . The pion is massless 
and the sigma particle has a mass of v^A/tt- The Goldstone theorem is satisfied. Lin and 



Serot |3T| have argued that the sigma meson should not be identified with the attractive 
s-wave interaction in the ir — n interaction, which is responsible for nuclear attraction. 
Rather, they argue that the sigma meson should have a mass which is at least 1 GeV if 
not more. This means that A is on the order of 50 or greater. 

The simplest approximation at finite temperature is the mean field approximation. 
One allows for v to be temperature dependent; hence the effective masses are temperature 
dependent as well. However, interactions among the particles or collective excitations are 
neglected. The pressure includes only the contribution of the condensate and of the 
thermal motion of the independently moving particles. Thus 

~ r In Z = ~ (fl -v 2 f + P (T, m a ) + (N - 1)P Q (T, m n ) . (9) 



V 4 

The pressure of a free relativistic boson gas can be written several ways: 

This is a relatively simple but surprisingly powerful first approximation which allows 
one to gain much insight into the behavior of relativistic quantum field theories at high 
temperature. It was used in all the pioneering papers. 
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One expects that as the temperature is raised, thermal fluctuations will tend to disor- 
der the condensate field v, and at sufficiently high temperature it may even disappear. If 
there is a second order phase transition then the correlation length should go to infinity, 
which is equivalent to the effective sigma mass going to zero. With such an expectation 
one may expand the free boson gas pressure about zero mass to obtain 

N 7T 2 . m 2 T 2 m 3 T , . 

Po(T, m) = — T + + ■•-. 11 

uv ' ; 90 24 12tt v ; 

Since the masses are proportional to the square root of A it is generally inconsistent to 
retain the cubic term in m because there exist loop diagrams which are not included in 
the mean field approximation but which contribute to the same order in A. Therefore we 
take 



P(T, v) = N—T A + —v 2 
v ' ; 90 2 



2 _ AT + 2 - 
J* 12 



~\v\ (12) 



where the pion and sigma masses have been expressed in terms of A, v and f n . Maximizing 
the pressure with respect to v gives 

v 2 = fl~^T\ (13) 

This result is easily understood. Going back to Eq. fl7|) we can differentiate InZ with 
respect to v with the result that 

v 2 = fl ~ 3<a 2 > - (vr 2 ) (14) 
as long as we choose (a) = 0. For any free bosonic field <fr with mass m 

M=f$r>h+h (15) 

where u = \/p' 2 + m 2 . In the limit that the temperature is greater than the mass (0 2 ) — > 
T 2 /12. This yields directly Eq. ©. 

The condensate goes to zero at a critical temperature given by 

T 2 = ^—f. (16) 

Above this temperature thermal fluctuations are too large to allow a nonzero condensate. 
It is a straightforward exercise to show that the pressure and its first derivative is con- 
tinuous at T c but that the second derivative is discontinuous. This is therefore a second 
order phase transition. 

There are two major problems with the mean field approximation as described. The 
first is that the pion has a negative mass squared at every temperature greater than zero. 
Not only is the Goldstone theorem not satisfied, but there are tachyons as well! The 
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sigma particle also gets a negative mass squared at temperatures above ^J8/(N + 2) f n < 
T c . This violation of basic physical principles is resolved by recognizing that the finite 
temperature corrections to the squared masses are proportional to AT 2 , and that one loop 
self-energy corrections, not included in the mean field analysis, are of the same order. This 
can be understood with the following analysis. 

At high temperatures, when the masses can be neglected in the loops, the mean field 
result is obtained by combining Eqs. @ and (p"3|). 

2 N + 2 ^2 

ml = 2A/ 2 -^AT 2 (17) 

The full one loop self-energies for pions and the sigma meson are drawn in Figures 1 and 
2. If one chooses (a) = then there are no 1PR diagrams and the tadpoles should not 
be included; they are already included in the temperature dependence of v. One may 
check this by fixing v = f n and then computing the tadpole contributions to the effective 
masses. One gets precisely Eq. ([T?D . The diagrams involving the 4-point vertices con- 
tribute an amount (N + 2)AT 2 /12 to both the pion and sigma meson self-energies. When 
evaluated in the high temperature approximation and at low frequency and momentum 
the 1PI diagrams involving the 3-point vertices may be neglected. (This follows from 
power counting. These diagrams involve two propagators instead of one, and so are only 
logarithmically divergent in the UV in the vacuum. The other diagrams are quadratically 
divergent, which leads to a T 2 behavior at finite temperature.) When all contributions of 
order AT 2 are included the pole positions of the pion and sigma propagators move with 
the result that below T c 

ml = 774 + 11^ = 

m 2 +n CT = 2A/ 2 (l-T 2 /T 2 ) , (18) 

-A/. 2 + II* = ^A (T 2 - T 2 ) . (19) 

The Goldstone Theorem is satisfied, there are no tachyons, and restoration of the full 
symmetry of the Lagrangian above T c is evident. 

The second major problem is that long wavelength fluctuations very near the phase 
transition cannot be treated with perturbation theory because the self-interacting boson 
fields become massless just at the transition. Although this is a well-known problem 
in the statistical mechanics of second order phase transitions, exactly how it affects the 
critical temperature is not known for the linear sigma model in 3 + 1 dimensions. This 
is a topic for further study. The result presented here must be accepted for what it is: a 
one loop estimate of the critical temperature. 



and above T c 



ml 



222 
m„ = m„ = m* 
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2.2 Nonlinear sigma model 

The nonlinear sigma model may be defined by the Lagrangian 

C = \ (W 2 (20) 

together with the constraint 

fl = &{^t). (21) 

The partition function is 

Z = J 5 [fl - $ 2 ) exp j jf" dr jf d 3 a; £ j . (22) 

Because the length of the chiral field is fixed and cannot be changed by thermal fluctua- 
tions it is often said that chiral symmetry breaking is built into this model and therefore 
there can be no chiral symmetry restoring phase transition. On the other hand, the linear 
sigma model does undergo a symmetry restoring phase transition. Taking the quartic 
coupling constant A to infinity essentially constrains the length of the chiral field to be 
f n just as in the nonlinear model. The critical temperature, however, is independent of 
A at least in the mean field approximation. So it would seem that the phase transition 
survives. If this is true, then one ought to be able to derive it entirely within the context of 
the nonlinear model. That is what we shall do, although it involves a lot more effort than 
treatment of the linear model in the mean field approximation. Since the only parameter 
in the model is /„-, and we are interested in temperatures comparable to it, we cannot 
do an expansion in powers of T/ f n . The only other parameter is N, the number of field 
components. This suggests an expansion in 1/N. 

Begin by representing the field-constraining delta function by an integral. 



J [d*] [db'\ exp IfdrJ d 3 x [C + ib' (V - /*)] 1 (23) 



As with the linear model, we define the first N — 1 components of $ to be the pion field 
and the last component to be the sigma field. We allow for a zero frequency and zero 
momentum condensate of the sigma field referred to as v. Following Polyakov we 
also separate out explicitly the zero frequency and zero momentum mode of the auxiliary 
field V ' . Integrating over all the other modes will give us an effective action involving the 
constant part of the fields. We will then minimize the free energy with respect to these 
constant parts, which is a saddle point approximation. Integrating over fluctuations about 
the saddle point is a finite volume correction and of no consequence in the thermodynamic 
limit. The Fourier expansions are 

$,(x,r) = yr^rH ££e^+^ >7r 4 (p,n), 
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$jv(x,r) = v + <j(x,t) =v + 



^(X'P+WnT) 



a(p,n) , 



i(x'P+j/„r) 



b(p,n) . 



(24) 



One must remember to exclude the zero frequency and zero momentum mode from the 
summations. The field $ must be periodic in imaginary time for the usual reasons, but 
there is no such requirement on the auxiliary field b, hence u n = lunT and v n = imT . 
Since the field b has dimensions of inverse length squared we inserted another factor of T 
so as to make its Fourier amplitude dimensionless, as they are for the other fields. The 
action then becomes 



S = [drj^x^- [(^) 2 - m V + (^f- m V] - tb (2va + tt 2 + a 2 )} 



+ - m < 



(25) 



Note that terms linear in the fields integrate to zero because (7Tj) = (a) = (b) = 0. 

An effective action is derived by expanding exp(S') in powers of b and integrating over 
the pion and sigma fields. The term linear in b vanishes on account of 6(0, 0) oc (b) = 0. 
The term proportional to b 2 is not zero and is exponentiated, thus summing a whole 
series of contributions. The term proportional to b 3 is not zero either and it, too, may be 
exponentiated, summing an infinite series of higher order terms left out of the order b 2 
exponentiation. After making the scaling b — > bj \JlN the effective action becomes 

S'eff = -- ( w « + P 2 + m2 ) [ # (P> n ) ■ *(-P> ~ n ) + ^(P> n)a{-p, -n)} 



n p 



n(p, uj n ,T,m) + — 



N u 2 n + p 2 + m' : 



+ -m' 



{fl - v 2 ) pv + o (lyVN) 



Note that only even Matsubara frequencies contribute in the 6-field: v, 
may have been anticipated. There appears the one loop function 



6(p,2n)6(-p,-2n) 

(26) 

2t™T. This 



n(p, u n , T, m) = T Y J 



d 3 k 



(2tt) 3 (uj n - oji) 2 + (p - k) 2 + m 2 lu 2 + k 2 + m 2 



(27) 



The effective action is an infinite series in b. The coefficients are frequency and momentum 
dependent, arising from one loop diagrams. The coefficient of the term quadratic in b in 
S c fi is illustrated schematically in Figure 3. In addition, each successive term is suppressed 
by 1/ VN compared to the previous one. This is the large expansion. 
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The propagators for the pion and sigma fields are of the usual form 



D (p,w n ,m) 



2,2, 2 
U n + V + m 



(28) 



with an effective mass m yet to be determined. The propagator for the 6-field is more 
complicated, being 



D b (p,u n ,m) = U(p,u n ,T,m) + 



N lu 2 + p 2 + m 2 



(29) 



The value of the condensate v is not yet determined either. 

Keeping only the terms up to order b 2 in S c g (the rest vanish in the limit iV — > oo) 
allows us to obtain an explicit expression for the partition function and the pressure. This 
includes the next to leading order in N. 



2 



- -TV 
2 ^ 



(2tt) 3 
d 3 p 



In 



II(p, u n ,T,m) + 



N u 2 + p 2 + m 2 



(30) 



The second term under the last logarithm should and will be set to zero at this order. It 
may be needed at higher order in the large N expansion to regulate infrared divergences. 

The pressure is extremized with respect to the mass parameter m. Therefore dPj dm 2 = 
0. From the initial expression for Z this is seen to be equivalent to the thermal average 
of the constraint. 

fl = <* 2 > = v 2 + (* 2 ) + <^ 2 > (31) 

If an approximation to the exact partition function is made, like the large iV expansion, 
this constraint should still be satisfied. It may, in fact, single out a preferred value of m. 

To leading order in iV we may neglect the term involving II entirely. The pressure is 
then 

P = im 2 (ft - v 2 ) +NP (T,m) . (32) 

The pressure must be a maximum with respect to variations in the condensate v. This 
means that 

dP/dv = -m 2 v = , (33) 
which is equivalent to the condition that (a) = 0. There are two possibilities. 

I. m = 0: There exist massless particles, or Goldstone bosons, and the value of the con- 
densate is determined by the thermally averaged constraint. This is the symmetry- 
broken phase. 
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2. v — 0: The thermally averaged constraint is satisfied by a nonzero temperature- 
dependent mass. There are no Goldstone bosons. This is the symmetry-restored 
phase. 

Evidently there is a chiral symmetry restoring phase transition! 

In the leading order of the large N approximation the particles are represented by 
free fields with a potentially temperature-dependent mass m. For any free bosonic field 



d 3 p 1 



dP (T, m)/dm 2 = m = I ^ - ^—^ (34) 



with u = y/p 2 + m 2 . Thus extremizing the pressure with respect to m is equivalent to 
satisfying the thermally averaged constraint. 

fl =v 2 + (tt 2 ) + (a 2 ) (35) 

Note however that the pion and sigma fields have the same mass and therefore (tz 2 ) = 
(N — l)(cr 2 ). Consider now the two different phases. 

In the asymmetric phase the mass is zero. The constraint is satisfied by a temperature- 
dependent condensate. 

NT 2 

v\T) = 11-^ (36) 
This condensate goes to zero at a critical temperature 

12 

T 2 = -jyf% (leading iV approximation) . (37) 

Exactly at T c the thermally averaged constraint is satisfied by the fluctuations of N 
massless degrees of freedom without the help of a condensate. 

In the symmetric phase the condensate is zero. The constraint is satisfied by thermal 
fluctuations alone. 

Thermal fluctuations decrease with increasing mass at fixed temperature. The constraint 
is only satisfied by massless excitations at one temperature, namely, T c . At temperatures 
T > T c the mass must be greater than zero. Near the critical temperature the mass should 
be small, and the flucutations may be expanded about m = as 



fl = NT 2 



1 m m 2 f m \ m 2 



1 2 AttT 8tt 2 T 2 ln 1 47T77 16n 2 T 2 + ' ' 



(39) 



As T approaches T c from above, the mass approaches zero like 

7T 

3T 



m(T) = -(T 2 -T 2 )+.... (40) 
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This is a second order phase transition since there is no possibility of metastable super- 
cooled or superheated states. 

The mass must grow faster than the temperature at very high temperatures in order 
to keep the field fluctuations fixed and equal to / 2 . Asymptotically the particles move 
nonrelativistically. This allows us to compute the fluctuations analytically. We get 



/ T \ 3/2 



-m/T 



(41) 



This is a transcendental equation for m(T). It can also be written as 



m 



T In 



NT mT 



Roughly, the solution behaves as 



m ~ T In (T 2 /T 2 ) 



(42) 



(43) 



It is rather amusing that, in leading order of the large N approximation, the elemen- 
tary excitations are massless below T c , become massive above T c , and at asymptotically 
high temperatures move nonrelativistically. 

The result to first order of the large N expansion provides good insight into the nature 
of the two-phase structure of the nonlinear sigma model, but it is not quite satisfactory 
for two reasons. First, it predicts N massless Goldstone bosons in the broken symmetry 
phase when in fact we know there ought to be only N — 1. Second, the square of the 
critical temperature is 12 f 2 /N whereas it is 12f 2 /(N + 2) in the linear sigma model in 
mean field approximation; we expect them to be the same in the limit A — > oo. Both these 
problems can be rectified by inclusion of the next to leading order term in N, namely, the 
contribution of the 6— field. 

It is natural to expect that the 6-field will contribute essentially one negative degree 
of freedom to the T 4 term in the pressure so as to give N — 1 Goldstone bosons in the 
low temperature phase. Therefore we move one of the N degrees of freedom and put it 
together with the b contribution as 



\m 2 {ft ~ v 2 ) - 
d 3 p 



N 



TJ2 



(2vr) £ 



In 



d 3 p 

(3 2 (uol+p 2 + m 2 )n 



In 



/3 2 (c 2 +p 2 + m 2 ; 



(44) 



The function U(p, cu n , T, m) can be reduced to a one dimensional integral 



n 



-/ 

T 2 V JO 



00 dkk 



8n 2 p 



In 



~k 2 +pk + A 2 ~ 
k 2 -pk + A 2 



(45) 
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where 



A 2 = A 2 (p,u n ,m) 



(46) 



but unfortunately cannot be simplified any further. In any case, to the order in iV to 
which we are working, the pressure is 

P = l -m 2 [fl - v 2 ) + (N - 1) P (T, to) + P(T, to) . (47) 

This can be thought of, in the low temperature phase, as iV — 1 Goldstone bosons with 
an interaction term P\. 

Because of the logarithm the main contribution to the interaction pressure will come 
when II is very small compared to one. This corresponds to very large values of the 
parameter A; in other words, to very high momentum, Matsubara frequency, or mass. In 
this limit 



n 



1 



4tt 2 A 2 



00 dkk 2 



h 3 (m/T) T 2 
A 2 



(48) 



uj ePv-l 4tt 2 

This may be thought of as a kind of high energy approximation, and we shall henceforth 
refer to it as such. Then 



d 3 p 



2 (2tt 



z n J 

with the dispersion relations 



(2tt) 3 
d 3 p 



In 



In 



/3 2 (u 2 n + P 2 + to 2 ) n" 

h 3 (u 2 n +p 2 )(u 2 n +p 2 + m 2 ) 

7T 2 (^ 2 +^ 2 )K 2 +^ 2 ) 



u±=p 2 + 2to 2 ± 2m\Jp 2 + to 2 . 
The interaction pressure can now be determined in the usual way to be 



(49) 



(50) 



Pi 



J 2tt 3 I 



-In 



(2vr)3 

1 _ e -/3^+(p) 



+ In 



-/Mp) 



-In 



l_ e -^-(p)]} 



(51) 



Note that h 3 (m/T) has no effect within this approximation. Note also that in the broken 
symmetry phase where to = the contribution of the 6-field cancels one of the massless 
degrees of freedom to give — 1 Goldstone bosons. 

Now we are prepared to examine the behavior of the system near the critical temper- 
ature with the inclusion of next-to- leading terms in N. We do an expansion in m/T as 
before. The pressure is, up to and including order to 3 : 



71 



(N - 1) — T 4 
V ; 90 



N + 2 
24 



-m 



'■T 2 + 1 -m 2 (f 2 ^v 2 ) + -^-m 3 T. 



(52) 
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In the high temperature phase where v — maximization with respect to m yields 



J TV 



N + 2 N m 



12 4ttT 



(53) 



This gives the same critical temperature as in the mean field treatment of the linear sigma 
model. 

12 

T c 2 = — — - fl (subleading N approximation) . (54) 
The mass approaches zero from above like 

MT) = =£±2 (I- - 7?) . (55) 

We leave it as an exercise for the reader to compute the asymptotic behavior of the mass 
with the inclusion of the subleading terms in N. 

The results obtained immediately above used an approximation for II which we re- 
ferred to as a high energy approximation. Relaxing this approximation can be done albeit 
at the cost of a numerical calculation. We do not attempt that in this paper. Of course, 
one should also go beyond the mean field approximation in the linear model. 



3 f n at Low Temperature 

Consideration of correlation functions at finite temperature is more involved than at 
zero temperature. Lorentz invariance is not manifest because there is a preferred frame 
of reference, the frame in which the matter is at rest. Thus spectral densities and other 
functions may depend on energy and momentum separately and not just on their invariant 
s. Also, the number of Lorentz tensors is greater because there is a new vector available, 
namely, the vector = (1,0,0,0) which specifies the rest frame of the matter. 

For a given four-momentum q it is useful to define two projection tensors. The first 
one, Pj? , is both three- and four-dimensionally transverse, 

p ij = gij _ <M_ ^ (56) 

with all other components zero. The second one, Pj^ , is only four-dimensionally trans- 
verse, 

pr = - (g^ - ^f- + P^j ■ (57) 

The notation is L for longitudinal and T for transverse with respect to q. There are no 
other symmetric second rank tensors which are four-dimensionally transverse. 
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In the usual fashion [RBI, 1331 one may construct a Green function for the axial vector 
current 

G£M = r —PZ(^<1), (58) 



where the spectral density tensor is 



fCM = ^E( 27r ) 3 ^-^ + ^)5(q-Pm + Pn) 

m,n 

x (e- s "/ T - e~ E ™ /T ) (n\A£(0) \m) (m\A u b (0) \n) . (59) 

The summation is over a complete set of energy eigenstates. The retarded, advanced and 
Matsubara Green functions are 

G^(9b,q) = + ie,q), (60) 

G^ v (q ,q) = G%(q -ie,q), (61) 
Gln^q) = G%(iu n ,q), (62) 

where e — ► + . 

Due to current conservation the spectral density tensor can be decomposed into lon- 
gitudinal and transverse pieces |35] . 

= [p L A^)Pr + p T A (q)Pr} (63) 

In general the spectral densities depend on q° and q separately as well as on the tempera- 
ture. In the vacuum we can always go to the rest frame of a massive particle, and in that 
frame there can be no difference between longitudinal and transverse polarizations, so 
that Pl = Pt = P- We also observe that P£ v + P^ v = —{g^ v — q fi q u /q 2 ). The pion, being 
a massless Goldstone boson, is special. It contributes to the longitudinal axial spectral 
density and not to the transverse one. In vacuum 

^fo) = [r~f ~ sr) p^ 2 ) + fttfW ■ (64) 

This may be taken to be the definition of the pion decay constant at zero temperature. 
In fact, one can write the pion's contribution as 

ffitfW = fh 2 W)PF ■ (65) 

This cannot be taken as the definition of the pion decay constant at finite temperature be- 
cause the contribution of the pion to the longitudinal spectral density cannot be assumed 
to be a 5 function in q 2 . In general the pion's dispersion relation will be more complicated 
and will develop a width at nonzero momentum. This smears out the delta function into 
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something like a relativistic Breit-Wigner distribution. Fortunately, the Goldstone Theo- 
rem |36| requires that there be a zero frequency excitation when the momentum is zero. 
(For a proof applicable to relativistic quantum field theories at finite temperature see |28| .) 
This implies that the width must go to zero at q = which results in a delta function 
at zero frequency. Explicit calculations support this assertion |37], |39| . Therefore it 
seems to make sense to define 

/^T)E21imr^p^ 0) q=0). (66) 

Physically this means that the pion decay constant at finite temperature measures the 
strength of the coupling of the Goldstone boson to the longitudinal part of the retarded 
axial vector response function in the limit of zero momentum. 

We shall study the pion's contribution to the spectral density only at temperatures 
small compared to f n . We shall study both the nonlinear and the linear sigma models. 
At low temperatures the sigma meson's contribution as a material degree of freedom is 
frozen out and one might expect the same dynamics to be operative in both models; in 
other words, one may expect the result to be the same and so independent of A. For tem- 
peratures approaching T c the problem is more difficult and is left for future investigation. 



3.1 Nonlinear sigma model 

The nonlinear sigma model was defined at the beginning of section 2.2. One can make 
a nonlinear redefinition of the field without changing the physical content of the theory. 
Various redefinitions may be found in the literature. We will first list the most common 
ones, and then we will compute fniT) for each of them, thereby illustrating that one 
always gets the same result. It is interesting to see how this comes about; it is also 
reassuring that it does. 

A convenient way to express the sigma and pion fields which explicitly contains the 
constraint is 



a 

TT 



U COS (0/ f w ) 

/^0sin(0// w ) 



(67) 



where <fi = \<f>\ and cf> = 4>l$- The Lagrangian may then be expressed in terms of the 
fields of choice. 

2 2 ^ 



1 

•j 



2 -<Vr • d»7T + 



1 (tt • d^iz) (tt • d^ir) 
2~ 



J TV 



71- 



sm 



J* 




(6£ 
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Another representation to consider is due to Weinberg M, who defines 



P = 




) 



7T 



(69) 




(70) 



1 df.p ■ d^p 

2 n+ v vm) 



(71) 



The (a, iz) representation is cumbersome because of the constraint, although it can 
be handled by the Lagrange multiplier method of section 2. However, it is inconvenient 
for exposing the physical particle content and for doing perturbation theory in terms of 
physical particles. Among the three physical representations we choose to work with here, 
it is interesting to note the range of allowed values of the fields. The magnitude of the 
p-field can range from zero to infinity, the magnitude of the 7r-field can range from 
to f n , and the magnitude of the (p-Held can range from to rrf^. This distinction is 
important when dealing with nonperturbative large amplitude motion; whether it makes 
any difference in low orders of perturbation theory is not known to us. 

The first step in our quest to extract the temperature dependence of f n from the 
theory is to obtain the form of the axial vector current in terms of the chosen fields. 
Starting from 



Every form of the axial vector current is an odd function of the pion field. 

Obviously it is not possible to compute the axial vector correlation function exactly. 
We will restrict our attention to low temperature. Roughly speaking, a loop expansion 
of the correlation function is an expansion in powers of T 2 / f%, with each additional loop 
contributing one more such factor. To one loop order we need the axial vector current to 



= -a d^ir + 7T <9 M cr 



(72) 



one directly computes 




(73) 
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third order in the pion field. 

TT 2 1 

A* = ^m 71 " + 97- fy w ~ t 77 I 71 " • fy 71 ") 
3o 2 1 

= -A^p+^^p-— p(p-^p) (74) 

We will also need the Lagrangian to fourth order in the pion field. 

1 

' r (0.^</>)(</>-^</>)-0 2 ^0-^0] 



A = ^ (TT • d,7r) (TT • 



6/ 2 



p 2 ^p-^p (75) 



4/ 2 

The correlation function (A^(x) Al(y)) will have a zero loop contribution from the pi- 
pi correlation function (d^ir 1 (x) d^rf (y)) , a one loop self-energy correction to the same 
pi-pi correlation function, and a one loop contribution from the correlation function 
{d^ l (x) ir j (y)ir k (y)dvir l (y)) involving four pions. These three contributions are illustrated 
in Figure 4. 

The contribution of the bare pion propagator D to the longitudinal spectral density 
is easily found to be 

p L A(qo,cO = f 2 w q 2 5(q 2 ) . (76) 

At zero temperature this is just the definition of the pion decay constant. 

The one loop pion self-energy may be computed by standard diagrammatic or func- 
tional integral techniques. The results are: 

T 2 

IM?) = —^pi 2 
n P (g) = (^-i)^ 2 

U^q) = i n7r (g) + jjll p (g). (77) 

These are quite dependent on the definition of the pion field! Nevertheless, it is worth 
noting that the Goldstone Theorem is satisfied on account of the fact that the self-energy 
is always proportional to q 2 . 

The final contribution comes from the correlation function of a pion at point x with 
three pions at point y. Again, standard diagrammatic or functional integral techniques 
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may be used. To express the answers, we gather together the contributions from the bare 
propagator, from the one loop self-energy, and from this correlation function, and quote 
the coefficient of the term f 2 q 2 5 (q 2 ) in the longitudinal part of the axial vector spectral 
density. 



7T 



p 







1 - 



T 2 



1 + (N- 
1 + {N- 



-(N-3) 



T 2 



12/. 2 



T 2 



24/ 2 
T 2 
'36/f 



8ft 



(N 



(78) 



In all three cases the results are the same and amount to a temperature dependence of 

N -IT 2 



i - 



12 ft 



It agrees with Eq. (|2|) for the only case which they can be compared: N 2 = 
calculation of Gasser and Leutwyler was verified by Eletsky and Kogan [40 . 



(79) 
N = 4. The 



3.2 Linear sigma model 

It is now not surprising to discover that the linear sigma model gives the same result for 
fn(T) at low temperature as the nonlinear sigma model. This is because the sigma meson 
is very heavy at low temperature and cannot contribute materially the way the pions do. 
However, the way in which it works out is very different. 

Let us go back to the axial vector current before shifting the sigma field. 

A,, = -a d^-K + 7T d^a (80) 

After making the shift a — ► v + a the current takes the form 

Aft = -v (9 M 7r - a d M 7r + 7r d^o . (81) 

By maximizing the pressure (minimizing the effective potential) with respect to v at each 
temperature we effectively sum all tadpole diagrams, leaving only 1PI diagrams in any 
subsequent perturbative treatment. If this is done, one's inclination is to identify v(T) 
with f n {T). This is wrong; f-^iT) has additional contributions, as we shall now see. 

The first contribution to f 2 (T) does come from v 2 (T) since it involves the cross term 
of <9 M 7r a (a:) with d u 7T a (y). Following the analysis of section 2.1, but at low temperature 
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rather than high, we simply leave out the contribution of the heavy sigma meson. This 
gives 



^ ~4 , A 2 



P(T,v) = (N-l)-T* + -v 



f2 N-l 



A 4 



Maximizing with respect to v gives 



v\T) = fl - 



N-l , 
12 



(82) 



(83) 



The T 2 / fl correction is identically the tadpole contribution to the vertex shown in Figure 
5. 

There is another, nonlocal, contribution to the vertex shown in Figure 5, corresponding 
to the emission and absorption of a virtual sigma meson. One might think that it is 
suppressed by the large sigma mass, m 2 = 2A/ 2 , but in fact this is compensated by the 
coupling constant A in the extra vertex. Evaluation of this diagram gives a contribution 
to/ 2 (T)ofT 2 /6/ 2 . 

Finally there is a contribution coming from the dressed pion propagator analogous to 
the nonlinear sigma model. The full one loop 1PI pion self-energy diagrams were already 
shown in Figure 1. We know that the sum of the momentum independent pieces is zero 
on account of Goldstone's Theorem. We only need the contribution which is quadratic in 
the energy and momentum of the pion. This can only arise from the so-called exchange 
diagram involving two aim vertices, also shown in Figure 5. In imaginary time (Euclidean 
space) it is 



n cx K, q ) = -4A 2 / 2 t]T I 



d 3 k 1 1 
(2^)3 u? + k 2 fa + uu n ) 2 + (k + q) 2 + ml 



(84) 



Since T <C m a it is easy to extract the piece quadratic in the momentum. Analytically 
continuing to Minkowski space (uj n — > iq ) it is g 2 T 2 /12/ 2 . 

The residue of the pion pole in the axial vector correlation function can now be 
obtained by adding the vacuum contribution, the pion self-energy correction, and the 
tadpole and nonlocal vertex corrections as follows. 



1 - 



N -IT 2 



+ 



IT 2 



The final result, 



fl{T) = fl 



12 fl Qfl 



N-2 T 2 
~^2~7? 



(85) 



is identical to that of the nonlinear sigma model. We remark that this cannot be used to 
compute the critical temperature since it was obtained under the condition that T <^ f n . 
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4 Scalar Condensate at Low Temperature 



The scalar condensate is defined as | (<&) |. Our convention has been to allow the last, iV'th 
component of the field to condense, and to refer to this as either v (if the field is shifted) 
or (a) (if the field is not shifted). In this section we use the latter convention. 

It is interesting to ask what happens to this condensate as a function of temperature 
in the nonlinear model. The constraint as an operator equation is / 2 = $ 2 and as a 
thermal average is / 2 = ($ 2 ); it is not f n = |(«&)|. The condensate indeed can change 
with temperature. In fact we can quite easily compute it to two loop order. Before doing 
so, we first discuss the connection of this condensate with the quark condensate (qq). 

In two flavor QCD one oftentimes associates the sigma and pion fields with certain 
bilinears of the quark fields. 

qq ~ a 



iqj 5 rq 



7T 



This association is made because the quark bilinears transform in the same way under 
SU(2)xSU(2) as the corresponding meson fields. The dimensions don't match so there 
must be some dimensionful coefficient relating them; this coefficient could even be a func- 
tion of the group invariant a 2 + 7r 2 ~ {qq) 2 — (qibTq) 2 . Does this particular combination 
of four-quark condensates change with temperature? The temperature dependence of the 
four-quark condensates at low temperatures was first calculated in f|l| with the help of 
the fluctuation-dissipation theorem. The contribution of pions alone was later discussed 
in j|2) using soft pion techniques. ^From |4l], f|2| one can read off the two condensates 



separately: 



and 



(iQ^rq) 2 ) 



1 - 



(W) 2 |0)-^(0|(g7 5 rg) 2 |0) 



12£ 



r- 



W 2 



T 2 



(0| (975-rgr |0) - — (0\(qqY\0). 



(86) 



(87) 



Therefore there is no correction to this group invariant to order T 2 // 2 inclusive. 



((qq) - (qisrqy) = (o| (my - (qi^qy |o> (88) 

This result is consistent with our analysis of the nonlinear sigma model in sections 2.2 
and 3.1. 

Now let us return to the business of computing the temperature dependence of the 
scalar condensate to one and two loop order. In terms of the three representations used 
in section 3.1 the sigma field is 
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(89) 



To second order in the pion field all three representations are the same. Using the free 
field expression for the thermal average of the field squared we get 



N — 1 / T 2 



.12/?. 



+ 



(90) 



For iV = 4, the only value for which we can quantitatively compare with QCD, this agrees 
with the result of Gasser and Leutwyler as quoted in Eq. @; it was also derived in an 
independent way by Eletsky p2|j . 

The coefficient of the term which is fourth order in the pion field differs in sign 
and magnitude among the three representations. It would be a miracle if the thermal 
average of \Jl — it 2 / ~f%, cos(0// 7r ) and the Weinberg expression were all the same! But 
regarding the order (T 2 /12/ 2 ) 2 we must recognize that the term which is second order in 
the pion field gets modified due to a one loop self-energy. This was computed for each 
representation in section 3.1 and the results listed in Eq. ([771). The term fourth order in 
the pion field can be evaluated using free fields. The result is 



<(0 2 ) 2 > = (iv 2 -i 




(91) 



and is obviously representation independent. The contributions for each representation 
are 




(92) 



where the second term in each line comes from the square of the pion field and the last 
term comes from the pion field in fourth order. The sum of all terms is identical in all 
three representations. 



(a)/U=l-(N-l) 



' T 2 



(N-l)(N-3) ( T 2 \ 2 
2 { 24/f J 



+ ■ ■ 



(93) 
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The miracle happens. It is a consequence of the fact that physical quantities must be 
independent of field redefinition. What is more, for N = 4 it agrees with the result of 
Gasser and Leutwyler quoted in Eq. (Q). However, we emphasize once more that this 
expression should not be used to infer a critical temperature because it was derived under 
the assumption that the temperature is small compared to f n . 

5 Summary and Conclusion 

In this paper we have focussed on the linear and nonlinear versions of the sigma model 
based on the group 0(N) at finite temperature. Models of this kind are prototypes for 
physical theories, such as QCD and electroweak theory. Our main goal was to understand, 
both conceptually and mathematically, whether the nonlinear model has a symmetry 
restoring phase transition analogous to that of the linear model. We did show that the 
nonlinear model has a second order phase transition by making use of the 1/N expansion. 
To leading and subleading orders the critical temperature is even the same as in the linear 
model. (This cannot be true in general; eventually there must be some dependence in 
the linear model on the value of the quartic coupling A.) This expansion was facilitated 
by the introduction of a Lagrange multiplier field. In this way we could see that there 
is a condensate at low temperature; this condensate decreases in just the right way so 
as to conserve the constraint on the field vector. There is one particular temperature 
for which the thermally averaged constraint is satisfied with no condensate and with all 
excitations massless. This is the critical temperature. We had to make a mathematical 
approximation at the subleading order to get an analytical result. We referred to this as a 
"high energy approximation" . It is directly analogous to what one does in the mean field 
approximation to the linear model. It would be interesting to relax this approximation; 
this is left as a future project. 

Another goal was to carefully define and show how to compute the "pion decay con- 
stant" and the "scalar quark condensate" at finite temperature within the scope of these 
models. The definitions also apply to full QCD but, of course, the results will generally 
be different. Only at very low temperatures and for N = 4 will the results be directly 
applicable to QCD for the reasons discussed in the introduction. Even within the context 
of the sigma models, however, it would be interesting to compute the next order correc- 
tion to f%(T) at low temperature. It would also be interesting to compute f%(T) near T c . 
These computations are now underway. 

In this paper we have not computed anything more complicated than a one loop 
diagram. Even the calculation of the "scalar quark condensate" to order (T 2 /12/ 2 ) 2 only 
required knowledge of one loop diagrams. What other interesting physical quantities can 
the reader compute in these models to one loop order? 

Natural extensions of these models to better approximate full QCD may be envisioned. 
Following the philospophy of chiral perturbation theory one may include higher derivative 
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terms in the Lagrangian. One may also add other mesonic and baryonic fields, especially 
the vector mesons. However, no matter how many extra terms are added, one is still 
restricted from discussing the quark-gluon plasma. 

If the Higgs particle turns out to have an exceptionally large mass, then a reasonable 
first approximation to the electroweak phase transition might begin with a nonlinear 
version of the Glashow-Weinberg-Salam model. Gauging the nonlinear sigma model 
would be a step in this direction. This topic is also under investigation. 

We hope to have stimulated the reader to make further progress on these very inter- 
esting topics at finite temperature! 

Acknowledgements 

We are grateful to A. Kovner and S. Jeon for discussions. This work was supported by 
the U.S. Department of Energy under grant number DE-FG02-87ER40328. 



26 



References 

[1] K. G. Wilson, Phys. Rev. D 7, 2911 (1973). 

[2] S.-K. Ma, Modern Theory of Critical Phenomena (Benjamin/Cummings, 1976). 

[3] D. J. Amit, Field Theory, the Renormalization Group, and Critical Phenomena (2nd 
ed., World Scientific, Singapore, 1984). 

[4] M. Gell-Mann and M. Levy, Nuovo Cimento 16, 705 (1960). 

[5] J. F. Donoghue, E. Golowich and B. R. Holstein, Dynamics of the Standard Model 
(Cambridge University Press, Cambridge, England, 1992). 

[6] S. Weinberg, Phys. Rev. 166, 1568 (1968). 

[7] S. Weinberg, Physica 96A, 327 (1979). 

[8] R. D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984). 

[9] F. Wilczek, Int. J. Mod. Phys. A7, 3911 (1992). 

[10] R. D. Pisarski, Phys. Lett. B 110, 155 (1982). 

[11] K. Rajagopal and F. Wilczek, Nucl. Phys. B399. 395 (1993); ibid. B404, 577 (1993). 

[12] Electroweak Physics and the Early Universe, NATO ASI Series B: Physics Vol. 338, 
ed. J. C. Romao and F. Freire (Plenum, New York, 1994). 

[13] D. A. Kirzhnits and A. D. Linde, Phys. Lett. B 42, 471 (1972). 

[14] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

[15] S. Weinberg, Phys. Rev. D 9, 3357 (1974). 

[16] For updates on the latest lattice results at finite temperature see the proceedings 
of the Lattice Conferences, the most recent being Nucl. Phys. B (Proc. Suppl.) 34 
(1993) and 42 (1995), or the proceedings of the Quark Matter Conferences, the most 
recent being Nucl. Phys. A 566 (1994) and A 590 (1995). 

[17] G. Boyd, S. Gupta, F. Karsch and E. Laermann, Z. Phys. C 64, 331 (1994). 

[18] E. Shuryak, Comm. Nucl. Part. Phys. 21, 235 (1994). 

[19] H. Leutwyler, in QCD 20 Years Later, Vol. 2, ed. P. M. Zerwas and H. A. Kastrup 
(World Scientific, Singapore, 1993). 

[20] F. Karsch, Phys. Rev. D 49, 3791 (1994). 



27 



[21] A. Kocic and J. Kogut, Phys. Rev. Lett. (1995). 

[22] M. Dine and W. Fischler, Phys. Lett. B 105, 207 (1981). 

[23] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987). 

[24] P. Gerber and H. Leutwyler, Nucl. Phys. B321, 387 (1989). 

[25] R. Dashen, S. Ma and H. J. Bernstein, Phys. Rev. 187, 345 (1969). 

[26] J. I. Kapusta, in Trends in Theoretical Physics, Vol. 2, ed. P. J. Ellis and Y. C. Tang 
(Addison Wesley, 1991). 

[27] G. Baym and G. Grinstein, Phys. Rev. D 15, 2897 (1977). 

[28] J. I. Kapusta, Finite Temperature Field Theory (Cambridge University Press, Cam- 
bridge, England, 1989). 

[29] V. F. Miiller and W. Riihl, Annals of Physics 168, 425 (1986). 

[30] P. Dreher, Phys. Lett. B 281, 127 (1992). 

[31] W. Lin and B. D. Serot, Nucl. Phys. A512, 637 (1990). 

[32] A. M. Polyakov, Gauge Fields and Strings (Harwood, Chur, Switzerland, 1987). 
[33] E. S. Fradkin, Proc. Lebedev Physics Institute 29, 6 (1965). 

[34] A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, 
McGraw-Hill (1971). 

[35] J. I. Kapusta and E. V. Shuryak, Phys. Rev. D 49, 4694 (1994). 

[36] J. Goldstone, Nuovo Cimento 19, 154 (1961). 

[37] J. Goity and H. Leutwyler, Phys. Lett. B 228, 517 (1989). 

[38] A. Schenk, Nucl. Phys. A363, 97 (1991); Phys. Rev. D 47, 5138 (1993). 

[39] C. Song, Phys. Rev. D 49, 1556 (1994). 

[40] V. L. Eletsky and I. I. Kogan, Phys. Rev. D 49, R3083 (1994). 

[41] A. Bochkarev and M. Shaposhnikov, Nucl. Phys. B268, 220 (1986). 

[42] V. L. Eletsky, Phys. Lett. B 299, 111 (1993). 



28 



Figure Captions 



Figure 1: One loop self-energy diagrams for the pion in the linear sigma model. The 
dashed lines represent pions and the solid lines represent sigma mesons. The three-point 
vertices are —Xv and the four-point vertices are —A/4. The = and 7^ indicate that the 
pion in the loop has the same or different quantum number than the external pion, re- 
spectively. If v is fixed at its vacuum value of f v then the two tadpoles contribute. If v is 
allowed to vary with temperature by maximizing the pressure then the tadpoles are not 
to be included in the self-energy; their effect is already included in the mean field mass 
via v(T). 

Figure 2: One loop self-energy diagrams for the sigma meson in the linear sigma model. 
See Figure 1 for remarks. 

Figure 3: Contribution to the effective action at finite temperature in the nonlinear sigma 
model corresponding to Eq. the wavy lines represent the Lagrange multiplier field 

b. 

Figure 4: Vertex and self-energy contributions to the axial vector correlation function in 
the nonlinear sigma model. 

Figure 5: Vertex and self-energy contributions to the axial vector correlation function in 
the linear sigma model. See Figure 1 for remarks. 
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